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Abstract Periodically modulated nonlinear oscillators often display bistability of 
forced vibrations. This bistability can be used for new types of quantum measure- 
ments. They are based on switching between coexisting vibrational states. Since 
switching is accompanied by a large change of the amplitude and phase of forced vi- 
brations, the measurements are highly sensitive. Quantum and classical noise plays 
dual role. It imposes a limitation on sensitivity in the familiar regime of a bifurca- 
tion amplifier. On the other hand, it makes it possible to use a bistable modulated 
oscillator in a new regime of a balanced dynamical bridge. We discuss the switching 
probabilities and show that they display scaling with control parameters. The criti- 
cal exponents are found for different types of bifurcations and for different types of 



1 Introduction 

Bistability of vibrational states in modulated systems and fluctuation-induced switch- 
ing between these states have attracted much attention recently. Experiments have 
been done on such diverse systems as electrons H) and atoms 12 O in modu- 
lated traps, rf-driven Josephson junctions JUS, and nano- and micromechanical 
resonators EH SHI. These systems have small vibration damping, the quality fac- 
tor may reach 10 4 — 10 5 . Therefore even a comparatively small resonant field can 
lead to coexistence of forced vibrations with different phases and amplitudes. This 
is illustrated in Fig.Q] which refers to a simple model relevant to many of the afore- 
mentioned experiments: an underdamped nonlinear classical oscillator driven close 
to resonance, with equation of motion 



Mark Dykman 

Department of Physics and Astronomy, Michigan State University, East Lansing, MI 48824, USA 
e-mail: c!ykman @pa .msu . edu| 



1 



2 



M. I. Dykman 



q + C0Qq + Yq 3 + 2rq=AcosC0Ft. (1) 

Here, q is the oscillator coordinate; Oo is its eigenfrequency, F is the friction coeffi- 
cient, r <C (Oo, and 7 is the nonlinearity parameter. The frequency of the modulating 
field (Of is assumed to be close to (Oq. In this case, for comparatively small modu- 
lation amplitude A, even where the oscillator becomes bistable its forced vibrations 
are nearly sinusoidal, q(t) = acos(©f? + 0) iflOl . 



Fig. 1 Bistability of a nonlin- 
ear Duffing oscillator {TJ. The 
solid lines show the squared 
amplitude of forced vibra- 
tions as a function of the 
squared amplitude of the driv- 
ing force. At the bifurcational 
values A B , 2 one of the stable 
vibrational states disappears. 



a 




The dynamic bistability of the oscillator is advantageous for measurements. The 
idea is to make the oscillator switch between the states depending on the value of 
the parameter to be measured. Switching leads to a strong change in the system that 
can be easily detected, leading to a high signal-to-noise ratio in a measurement. This 
has been successfully used for fast and sensitive measurements of the states of dif- 
ferent types of Josephson junction based qubits, including quantum non-demolition 
measurements ll4l [5l [TTI [T21l . 

So far the experiments were done in the bifurcation amplifier mode, where the 
control parameter is swept through a bifurcation point (for example, the field ampli- 
tude A was swept through Ab2, see Fig . |T]> . The position of the bifurcation point, i.e., 
the value of A where switching occurs, depends on the state of the measured qubit. 
However, instead of happening for a certain A, as expected from Fig. [T] switching 
happens at random within a certain parameter range of A near Ab2- 

The randomness of the switching field is a consequence of fluctuations in the 
oscillator. They lead to switching even before the control parameter reaches its bi- 
furcational value. This is analogous to activated switching out of a potential well 
studied by Kramers 1131 . However, in the case of an oscillator the stable states are 
not minima of a potential, and there is no static potential barrier that needs to be 
overcome. Switching of a modulated oscillator is an example of metastable decay 
of systems far from thermal equilibrium, the phenomenon of a broad interest. 

Theoretical analysis of metastable decay requires developing methods for cal- 
culating the decay probability and finding out whether decay displays any universal 
system-independent features, like scaling dependence on the control parameters. For 
classical systems, scaling of the decay rate was indeed found for systems close to a 
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bifurcation point, both in the cases of equilibrium Ifl4l [151 [T6l and nonequilibrium 
systems ifTTl [T8l [19| . In the latter case a scaling crossover may occur as the system 
goes from the underdamped to overdamped regime while approaching the bifurca- 
tion point [ 20J . Such crossover occurs also for quantum tunneling in equilibrium 
dissipative systems 1211 . 

In this paper we study decay of metastable vibrational states in quantum dissi- 
pative systems close to bifurcation points 1221 . This is necessary for understand- 
ing the operation of a modulated oscillator in the regime of a bifurcation amplifier. 
We show that at low temperatures decay occurs via quantum activation. This is a 
specific process that has no analog in systems in thermal equilibrium l23l l24l . As 
tunneling, quantum activation is due to quantum fluctuations, but as thermal activa- 
tion, it involves diffusion over an effective barrier separating the metastable state. 
As we show, near a bifurcation point quantum activation is more probable than tun- 
neling even for T — > 0. We find that the decay rate W scales with the distance to the 
bifurcation point 77 as | lnW| Jp. The scaling exponent is t, — 3/2 for resonant 
driving, cf. Eq. ([T). We also consider parametric resonance in a nonlinear oscilla- 
tor and show that in this case £, = 2. In addition, | InW | displays a characteristic 
temperature dependence. 



2 Quantum Kinetic Equation for a Resonantly Driven Oscillator 

The Hamiltonian of a resonantly driven nonlinear oscillator is 

Ho(t) = ^p 2 + ^q 2 +^Yq 4 -qAcos(co P t). (2) 

The notations are the same as in equation of motion £[), p is the oscillator momen- 
tum. We assume that the detuning 8(0 = (Of — <$o of the modulation frequency (Of 
from the oscillator eigenfrequency ©0 is small and that y8(0 > 0, which is necessary 
for bistability; for concreteness we set 7 > 0. 

It is convenient to switch from q,p to slowly varying operators Q,P, using a 
transformation 

q = Qes (Q COS COff+P sin (Opt), p = -C xes (0 F (gsin <X>Ft - PCOS (OFt ) (3) 

with C res = (8(Of8(o/3Y) 1 / 2 . The variables Q,P are the scaled coordinate and mo- 
mentum in the rotating frame. They are canonically conjugate, 

[p,Q] = -a, X = 3fry/8(0 F 8(0. (4) 

The parameter A plays the role of the effective Planck constant. We are interested in 
the semiclassical case; A is the small parameter of the theory, A <C 1. 

In the rotating wave approximation the Hamiltonian (O becomes Hq = (h/X) 8(0 g, 
with 



4 



M. I. Dykman 



§ = g (Q,P) = l(Q l + I*-l)*-p 1 / 2 Q, p = 3r A 2 /32(oj(8(o) 3 (5) 

(in the case y, 8(0 < one should use g— ► — g,Ho — > — (h/X)8cog). The function g is 
shown in Fig. [2] It plays the role of the oscillator Hamiltonian in dimensionless time 
T? = t\8co\. The eigenvalues of g give oscillator quasienergies. The parameter j3 in 
Eq. © is the scaled intensity of the driving field. For weak damping the oscillator is 
bistable provided < /3 < 4/27. The Heisenberg equation of motion for an arbitrary 
operatorM is M = dM/dx = —iX~ l [M,g]. 




Fig. 2 Quasienergy of a resonantly driven nonlinear oscillator g(Q,P) (left panel) and its cross- 
section by the plane P — (right panel). The plot refers to ji = 1 /270 in Eq. Thin horizontal 
lines in the right panel show (schematically) quasienergy levels for quantized motion around the 
local maximum of g{Q,P). In the presence of dissipation the states at the local maximum and 
the minimum of g(Q, P) become stable. They correspond, respectively, to forced vibrations with 
small and large amplitude a in Fig. [T] The arrows in the right panel show relaxation to the state 
of small-amplitude vibrations, tunneling from this state with constant quasienergy g, and quantum 
activation. The latter corresponds to quantum diffusion over quasienergy away from the metastable 
state, which accompanies relaxation I23II24I . 

We will consider two major relaxation mechanisms of the oscillator: damping 
due to coupling to a thermal bath and dephasing due to oscillator frequency mod- 
ulation by an external noise. Usually the most important damping mechanism is 
transitions between neighboring oscillator energy levels. They result from coupling 
linear in the oscillator coordinate. Since the energy transfer is « h(Oo, in the ro- 
tating frame the transitions look instantaneous. Phenomenologically, the resulting 
relaxation may be described by a friction force proportional to velocity, as in (Q~|). 
Microscopically, such description applies in the case of Ohmic dissipation, i.e., cou- 
pling to Ohmic bath. However, we do not have to assume that dissipation is Ohmic. 
The only assumption needed for the further analysis is that the density of states of 
the reservoir weighted with the interaction be smooth in the frequency range, which 
is centered at coq and has a width that largely exceeds F, |5t»|. 

We will assume that the correlation time of the noise that modulates the oscil- 
lator frequency is also short compared to l/|5o|, so that the noise is effectively 
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(5 -correlated in slow time T. Then the quantum kinetic equation is Markovian in the 
rotating frame. It has a familiar form (cf. l25l ) 

p = d T p = iX- l [p,g}-fp -fP h p, (6) 

where tp describes damping, 

fp =r\8co\~ 1 [(n+ l)(d f dp — 2ap&* + pd'a) +n(dd^ p -2a 1 pd + pdd^)] , (7) 

and r ph p describes dephasing, 

fp h P = n h \8co\- 1 [tfa, [fl + a,p]] . (8) 

Here, T and r ph are the damping and dephasing rates, d = (2X) l ^ 2 (Q + iP) is the 
oscillator lowering operator, and n = [exp (Ticoq / kT) — 1] _1 is the oscillator Planck 
number. In the classical mean-field limit one can obtain from (O, (0 the same equa- 
tion of motion as Eq. (Q]) in the rotating wave approximation. 
In what follows we use dimensionless parameters 

a = \Sto\/r, KP h = rP h /xr. (9) 

We assume that K-P h <l. This means that the intensity of phase fluctuations may be 
comparable to the intensity of quantum fluctuations associated with damping, which 
is oc XT, see below, but that T ph <C T. 

Metastable decay was studied earlier for additively and parametrically driven 
oscillators at T = K" ph = where there is detailed balance l26l l27l [28], and the 
lowest nonzero eigenvalue of p was studied numerically in Ref. (29) . However, the 
T = K ph = solution is fragile. It can change exponentially strongly already for 
extremely small T, JC ph 11231 l24l . The analysis l23l l24l revealed the mechanism of 
quantum activation over a quasienergy barrier, but the results referred to the case 
where the damping-induced broadening of quasienergy levels is small compared to 
the typical interlevel distance. This condition necessarily breaks sufficiently close to 
a bifurcation point where the level spacing becomes small as a consequence of the 
motion slowing down. 



3 Wigner Representation 

The analysis of metastable decay near a bifurcation point can be conveniently done 
in the Wigner representation, 

Pw(Q,P)= Jd^-^pfQ+ l -^Q- l -A, (10) 
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where p (2i : G2) = (Gi IP IG2} is the density matrix in the coordinate representation. 
Using Eqs. d4T>-(fTob one can formally write the equation for pw as a sum of terms 
proportional to different powers of the effective Planck constant A, 

pw = -V-(Kp w )+XL^p w + X 2 L^p w , (11) 

where K = (Kg,Kp) and V = (dg,dp). Vector K determines evolution of the density 
matrix in the absence of quantum and classical fluctuations, 

K Q = d P g-n- l Q Kp = -d Qg -Q- 1 P. (12) 

This evolution corresponds to classical motion 

Q = K Q , P = K p . (13) 

The condition K = gives the values of Q,P at the stationary states of the oscillator 
in the rotating frame. 

The term L* 1 ' in (fTTT i describes classical and quantum fluctuations due to damp- 
ing and dephasing, 



n+ 2 



V 2 + KV h (Qd P ~Pd Q ) 2 



(14) 



These fluctuations lead to diffusion in (Q,P)-space, as seen from the structure of 
LS 1 ': this operator is quadratic in <9g,<3/>. 

The term in (fTTT i describes quantum effects of motion of the isolated oscilla- 
tor, 

LW = ~(Qd P -Pd Q )V 2 . (15) 

In contrast to IS 1 ', the operator Z/ 2 ' contains third derivatives. Generally the term 
X 2 U 2 >p-w is not small, because pw varies on distances ~ X. However, it becomes 
small close to bifurcation points, as shown below. 



3.1 Vicinity of a Bifurcation Point 



From ( fT2b . ( fT3l ), for given reduced damping £2 1 the oscillator has two stable and 
one unstable stationary state in the rotating frame (periodic states of forced vibra- 
tions) in the range (.Q) < j3 < j8g 2 '(X2) and one stable state outside this range 
iflOl . with 

2 



3(1,2) 



27 



-9X2- 2 t(1-3X2- 2 ) 3/2 



(16) 



At jS^ 1 ' and /3^ 2 ' the stable states with large and small Q 2 
and small vibration amplitudes), merge with the saddle state (saddle-node bifurca- 



-P , respectively (large 
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tion). The values of Q,P at the bifurcation points 1, 2 are 



Q B = P B l/2 Y B (Y B -l), P B = p B 1/2 Q- 1 Y B (Y B = Ql+Pl), 



(17) 



with 



(1,2) 



2±(l-3f2- 2 ) 1 / 2 



(18) 



In the absence of fluctuations the dynamics of a classical system near a saddle- 
node bifurcation point is controlled by one slow variable |30| . In our case it can 
be found by expanding Kq p in 8Q = Q — Q B , 8P = P — P B , and the distance to 
the bifurcation point 77 = /3 — f5 B . The function Kp does not contain linear terms in 
8Q, 8P. Then, from 02), P slowly varies in time for small 8Q, 8P, 77. 

On the other hand 



-m-'iSQ-ciBSP), a B = Q(2Y B -\) 



(19) 



Therefore the relaxation time of Q is £2/2, it does not depend on the distance to the 
bifurcation point. As a consequence, Q follows P adiabatically, i.e., over time ~ £2 
it adjusts to the instantaneous value of P. 



4 Metastable Decay near a Bifurcation Point 

The adiabatic approximation can be applied also to fluctuating systems, and as we 
show it allows finding the rate of metastable decay. The approach is well known for 
classical systems described by the Fokker-Planck equation l3~ll . It can be extended 
to the quantum problem by factoring pw into a normalized Gaussian distribution 
over 8Q = 8Q — a B 8P and a function pw{8P) that describes the distribution over 
8P, 

(8Q-a B 8P) 2 



pw « const x exp 



-2 



A(2n+l)(l+a 2 



B) 



Pw(8P). 



In the spirit of the adiabatic approximation, p~w can be calculated disregarding small 
fluctuations of 8Q. Formally, one obtains an equation for pw by substituting the 
factorized distribution into the full kinetic equation (fTTT i and integrating over 8Q. 
This gives 

P~w ~ d P [PwdpU + A 2> B d P pw] , (20) 
where U and & have the form 

U= l -b(SPf-±l3 B l/2 i 1 SP, J] 



s 
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= £2 



-i 



(21) 



with 



1 -P 1 b /2 Q 2 (3Y b -2). 



In < T20b , d2TT i we kept only the lowest order terms in 8P,j5 — /3b, X. In particular we 
dropped the term — A 2 gfi<3pPw/4 which comes from the operator L 1 - 2 ^ in Eq. (fTTT i. 
One can show that, for typical \8P\ ~ |J7|'^ 2 , this term leads to corrections ~ T],A 
top w . 

Eq. d20b has a standard form of the equation for classical diffusion in a potential 
U(8P), with diffusion coefficient However, in the present case the diffusion 
is due to quantum processes and the diffusion coefficient is h for T — > 0. 



4.1 Scaling of the Rate ofMetastable Decay 

For r\b > the potential t/ d2Tb has a minimum and a maximum. They correspond 
to the stable and saddle states of the oscillator, respectively. The distribution pw has 
a diffusion-broadened peak at the stable state. Diffusion also leads to escape from 
the stable state, i.e., to metastable decay. The decay rate W is given by the Kramers 
theory O, 

R IX 2 1 /2| 7? |3/2 

W = CeT**! x , R A = -^ji , (22) 

with prefactor C = iz~ l (br\ /2) l l 2 fi B l ^\8(o\ (in unsealed time?)- 

The rate d22l displays activation dependence on the effective Planck constant 
X. The characteristic quantum activation energy Ra scales with the distance to the 
bifurcation point 77 = — /3g as rj 3 / 2 . This scaling is independent of temperature. 
However, the factor &b in Ra displays a characteristic T dependence. In the absence 
of dephasing we have = 1 /2X2 for « <C 1, whereas 2>b = kT /JicoqQ for « 3> 1. 
In the latter case the expression for W coincides with the result ifTTl . 

In the limit £2 ^> 1 the activation energy d22l) for the small-amplitude state has 
the same form as in the range of j3 still close but further away from the bifurcation 
point, where the distance between quasienergy levels largely exceeds their width 
ll23ll . We note that the rate of tunneling decay for this state is exponentially smaller. 
The tunneling is shown by the dashed line in the right panel of Fig. [2] The tunneling 
exponent for constant quasienergy scales as rj 5 / 4 |[T8l[32l . which is parametrically 
larger than 77 3 / 2 for small rj [for comparison, for a particle in a cubic potential OTT i 
the tunneling exponent in the strong-damping limit scales as T] l2D l. 

For the large-amplitude state the quantum activation energy ( l22l displays differ- 
ent scaling from that further away from the bifurcation point, where Ra°< j3 '^ 2 for 
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Q. 1 11231 . For this state we therefore expect a scaling crossover to occur with 
varying /3 . 



5 Parametrically Modulated Oscillator 

The approach to decay of vibrational states can be extended to a parametrically 
modulated oscillator that displays parametric resonance. The Hamiltonian of such 
an oscillator is 

Ho(t) = \p 2 + \q 2 H + Fcos((0 F 0] + l -yq A . (23) 

When the modulation frequency C0f is close to 2(flo the oscillator may have two 
stable states of vibrations at frequency (Of/ 2 (period-two states) shifted in phase by 
n iflOl . For F <C (Oq the oscillator dynamics is characterized by the dimensionless 
frequency detuning jj,, effective Planck constant A, and relaxation time £, 

cof(cof-2coq) 3\y\h F 

= f ' C ~2^r' (24) 

As before, A will be the small parameter of the theory. 

Parametric excitation requires that the modulation be sufficiently strong, £ > 1 . 
For such £ the bifurcation values of fi are 

4 1 ' 2) = T(1-C" 2 ) 1/2 , C>1- (25) 

If 7> 0, as we assume, for < [ijp the oscillator has one stable state; the amplitude 

of vibrations at g>f/2 is zero. As /i increases and reaches /ig this state becomes 
unstable and there emerge two stable period two states, which are close in phase 
space for small ji — fXg (a supercritical pitchfork bifurcation). These states remain 

(2) 

stable for larger /I. In addition, when /I reaches ji B the zero-amplitude state also 
becomes stable (a subcritical pitchfork bifurcation). The case y < is described by 
replacing jtx — > — fi. 

The classical fluctuation-free dynamics for fi close to Hb is controlled by one 
slow variable fl30ll . The analysis analogous to that for the resonant case shows that, 
in the Wigner representation, fluctuations are described by one-dimensional diffu- 
sion in a potential, which in the present case is quartic in the slow variable. The 
probability W of switching between the period-two states for small pL — /lg and the 

(2) 

decay probability of the zero-amplitude state for small /i — fig' have the form 



W = Cexp(-R A /X), fl A = |^|T7 2 /2(2«+l), 7? = jU - pi B (26) 
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(Hb = Us ' 2) )- The corresponding prefactors are = 2C ( B l) = 2 l l 2 K- x r C, 2 \pL B \ |jU - 
Hb\. Interestingly, dephasing does not affect the decay rate, to the lowest order in 
jj. — flB- This is a remarkable feature of quantum activation near bifurcation points 
at parametric resonance. 

From ( T26l l. at parametric resonance the quantum activation energy R A scales with 
the distance to the bifurcation point as r\ 2 . In the limit £ 3> 1 the same expression 
as d2oT i describes switching between period-two states still close but further away 
from the bifurcation point, where the distance between quasienergy levels largely 
exceeds their width. In contrast, the exponent for tunneling decay in this case scales 
as 7] 3 / 2 CD. 



6 Balanced Dynamical Bridge for Quantum Measurements 

Dynamic bistability of a resonantly driven oscillator can be used for quantum mea- 
surements in yet another way, which is based on a balanced bridge approach. As 
a consequence of interstate transitions there is ultimately formed a stationary dis- 
tribution over coexisting stable states of forced vibrations. The ratio of the state 
populations w\ ,w 2 is 

Wl /w 2 = W 2l /W l2 - exp[-(R A2 -R M )/X], (27) 

where W mn is the switching rate m — ► n. For most parameter values \R A \ — 1 ^ X, 
and then only one state is predominantly occupied. However, for a certain relation 
between j3 and £2, where R A \ ps R A2 , the populations of the two states become equal 
to each other. This is a kinetic phase transition 1331 . A number of unusual effects 
related to this transition have been observed in recent experiments for the case where 
fluctuations were dominated by classical noise ll9l [34l 1351 . 

In the regime of the kinetic phase transition the oscillator acts as a balanced 
dynamical bridge: the populations are almost equal with no perturbation, but any 
perturbation that imbalances the activation energies leads to a dramatic change of 
w\ t 2, making one of them practically equal to zero and the other to 1. Such a change 
can be easily detected using, for example, the same detection technique as in the 
bifurcation amplifier regime. An interesting application of a bridge is measurement 
of the statistics of white noise in quantum devices. The Gaussian part of the noise 
does not move a classical oscillator away from the kinetic phase transition. In con- 
trast, the non-Gaussian part, which is extremely important and extremely hard to 
measure, does, and therefore can be quantitatively characterized. 
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It follows from the above results that, both for resonant and parametric modulation, 
close to bifurcation points decay of metastable vibrational states occurs via quan- 
tum activation. The quantum activation energy Ra scales with the distance to the 
bifurcation point 77 as Ra « Jp, with | = 3/2 for resonant driving and ^ = 2 for 
parametric resonance. The activation energy Ra is smaller than the tunneling expo- 
nent. Near bifurcation points these quantities become parametrically different and 
scale as different powers of 77 , with the scaling exponent for tunneling (5 /4 and 3/2 
for resonant driving and parametric resonance, respectively) being always smaller 
than for quantum activation. 

The exponent of the decay rate Ra/X displays a characteristic dependence on 
temperature. In the absence of dephasing, for kT ^> ha>o we have standard thermal 
activation, R A °c 1 /T. The low-temperature limit is described by the same expression 
with kT replaced by h(Oo/2, 

Our results show that quantum activation is a characteristic quantum feature of 
metastable decay of vibrational states. Activated decay may not be eliminated by 
lowering temperature. It imposes a limit on the sensitivity of bifurcation amplifiers 
based on modulated Josephson oscillators used for quantum measurements ITTl[T2ll . 
At the same time, an advantageous feature of dynamically bistable detectors is that 
they can be conveniently controlled by changing the amplitude and frequency of the 
modulating signal. Therefore the value of the activation energy may be increased 
by adjusting these parameters. Another advantageous feature of these detectors, 
which has made them so attractive, is that they operate at high frequency and have 
a small response time, which can also be controlled. The results of this work apply 
to currently studied Josephson junctions, where quantum regime is within reach. 
They apply also to nano- and micromechanical resonators that can be used for high- 
sensitivity measurements in the regime of a balanced dynamical bridge. 
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